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Analysis of ray phase-space recirculation in an extended Budden problem 
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A three-wave Budden model with two resonance layers is constructed that allows recirculation of 
energy fluxes along a quadrangle in ray phase space. The transmission, reflection, and conversion 
coefficients for this extended Budden problem are calculated by ray phase-space methods and the 
modular-eikonal approach. Analytical and numerical results show that all coefficients exhibit inter- 
ference effects that depend on an interference phase calculated from the coupling constants and the 
area enclosed by the quadrangle. 
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The process of mode conversion, involving the interac- 
tion between two distinct waves propagating in a mag- 
netized plasma, represents an important mechanism by 
which heating and current drive by radio frequency may 
lead to a burning plasma. In a standard single-crossing 
mode-conversion scenario (see conversion region a in 
Fig. [1} , an incident ray from a primary (:r-propagating) 
wave propagates on its dispersion surface until it crosses 
transversely the dispersion surface of a secondary (k- 
propagating) wave. Mode conversion occurs when some 
of the energy flux associated with the primary wave is 
converted to energy flux for the secondary wave, which 
allows a secondary-wave ray to propagate away from the 
conversion region. The remaining energy flux associated 
with the primary wave is transmitted across the conver- 
sion region and allows a primary-wave ray to propagate 
away from the conversion region. The conservation of en- 
ergy (flux) dictates that the sum of the outgoing energy 
fluxes is equal to the incoming energy flux. 

The classic Budden problem is a double-conversion 
process in ray phase space [jj (see conversion regions b 
and c in Fig. [I]), whereby a primary incoming wave is 
converted to a secondary wave at a first conversion re- 
gion (b), which then propagates to a second conversion 
region (c) where it is converted to an outgoing (reflected) 
primary wave. Using ray phase-space methods and the 
modular-eikonal approach (2j, we extend the standard 
Budden problem by adding a second double-conversion 
process (see conversion regions a and d in Fig. [I]), and 
a second resonance layer, associated with a third wave 
possibly supported by an energetic-particle population 
in an inhomogeneous magnetized plasma. The extended 
three-wave Budden problem introduces the possibility of 
a leaky cavity being created in ray phase space that al- 
lows recirculation of energy fluxes. 

The modular-eikonal analysis of energy-flux recircula- 



tion around a leaky phase-space cavity in ray phase space 
was first performed by Brizard et al. 0, 0] , who consid- 
ered the linear interaction between two waves with dis- 
persion curves modeled by intersecting parabolas. In this 
earlier work, the modular-eikonal approach was success- 
fully applied to the double-crossing mode conversion pro- 
cess whereby a ray of one wave punctured the dispersion 
surface of another wave twice due to ray curvature (corre- 
sponding to the case of a near-tangential crossing). The 
ray phase-space analysis revealed the existence of four 
mode-conversion scenarios based on whether the wave en- 
ergies had equal or opposite signs, and whether the wave 
rays were co-propagating or counter-propagating in ray 
phase space (i.e., whether the signs of dk/dt were equal 
or opposite). While all four scenarios yielded transmis- 
sion and conversion coefficients exhibiting interference ef- 
fects associated with the area enclosed by the intersecting 
parabolas and the strength of the coupling constants at 
the two conversion points, only two scenarios involved 
recirculation of energy fluxes. 

The present work considers a modification of the three- 
wave Budden problem studied by Liang et al. 0, who 
considered the mode conversion in two-dimensional ray 
phase space (x, k) between an x-propagating wave with 
two positive-energy waves that are co-propagating along 
the fc-axis. The extended three-wave Budden problems, 
considered here and in Ref . 0] , introduce four conversion 
regions (a, 6, c, d) (instead of two for the classic Budden 
problem) that form the vertices of a quadrangle in ray 
phase space produced by the pair-wise intersections of 4 
uncoupled rays associated with three waves (see Fig. [TJ . 
The modular-eikonal analysis of Liang et al. showed 
that the reflection coefficient and one of the conversion 
coefficients exhibited interference effects that depended 
on the area enclosed by the quadrangle and the strength 
of the coupling constants at the four vertices. In addition, 
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FIG. 1: Counter-fc-propagating rays B± interacting with x- 
propagating ray A. Solid lines are dispersion curves repre- 
senting Eq. ([2]). The four conversion regions (a, b, c, d) form a 
quadrangle in ray phase space and arrows show the direction 
of wave-energy flow. The net propagation phase $ = <j> kdx 
is shown only in the boxed expressions. For each unit of en- 
ergy carried by the incoming left-propagating ray A (from 
x — + oo) introduced at conversion region a, the outgoing 
left-propagating ray A carries the transmitted energy flux T 
(left of conversion region b) , the outgoing up-propagating ray 
B+ carries the converted energy flux C+ (above conversion 
region c), the outgoing right-propagating ray A carries the 
reflected energy flux R (right of conversion region d) , and the 
outgoing down-propagating ray B- carries the converted en- 
ergy flux C_ (below conversion region a). Conservation of 
energy flux requires that T + C+ + C- + R = 1. 



Liang et al. [|J obtained excellent agreement between di- 
rect numerical integration of the coupled-wave equations 
and analytical results, thus confirming the validity of the 
modular-eikonal approach. 

The purpose of the present paper is to use the modular- 
eikonal approach to study the extended three-wave Bud- 
den problem in which two counter-propagating positive- 
energy waves are each interacting linearly with an x- 
propagating wave. This scenario allows us to investigate 
the validity of the modular-eikonal approach when energy 
fluxes are allowed to recirculate around a quadrangle in 
ray phase space. 

We begin with the extended three- wave Budden model 
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where the x-propagating wave A is represented by its 
dispersion function D A {k) = k 2 — k 2 (e.g., a magne- 
tosonic wave) and the two fc-propagating waves B± are 
represented by their dispersion functions D±(x,k) = 



(a±x)—j k. The bulk-ion and energetic-particle positive- 
energy waves are represented by the B + and £?_ waves, 
respectively, and the standard Budden problem is recov- 
ered with r;_ = 0. In Eq. ([I}, k = —id/dx and the 
strengths of the mode conversions A <-> B± are governed 
by the coupling constants rj±. The upward and down- 
ward fc-propagating waves (separated by the distance 2 a 
on the x-axis) are allowed to propagate (in x-space) by 
setting the parameter 7 > 0. By using the 3x3 dispersion 
matrix D(x, k) of Eq. |T]), the dispersion curves (shown 
in Fig. [T]) are obtained from the dispersion relation 
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0, (2) 



and the uncoupled dispersion curves cross in ray phase 
space at the conversion points (£( a ,d), fya,d)) = (c ± 
iko, Tko) and (x^ b} ,k^ b} ) = (-cr ±jk , ±k ). The 
energy-flux propagation is represented by the Hamil- 
ton equations (x = — dkD, k = d x D), which yield 
(x A , k A ) = (±2, 0) and (x±, k±) = (7, ±1), indi- 
cated by arrows in Fig. [TJ We note that the symmetric- 
quadrangle dispersion curves associated with the ex- 
tended Budden model (fTJ) facilitate the calculations done 
here, while the use of an asymmetric-quadrangle model 
(based on generic dispersion functions D± with X- 7^ x + 
and \k—\ 7^ \k+\) would not qualititatively change our 
results (as will be discussed below). 

The coupled differential equations of the extended 
Budden model Q: 



d 2 A(x) 
dx 2 



+ k 2 A(x) 
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dB±(x) , . „ . . 
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satisfy the wave energy-flux conservation law 



- ( j A (x) + x: ) = 0, 



(4) 



3=± 



where J A {x) = lm(A* dA/dx) and J±(x) = \l\B±\ 2 . 
The choice of the x-space propagation parameter 7 for 
waves B± does not affect the results of the modular- 
eikonal approach Q presented for our symmetric- 
quadrangle model. When 7 = 0, we replace B± with 
— [r)±/((T ± a;)] A in the wave equation for A{x) and we 
obtain the extended Budden equation 



d 2 A 
dx 2 



■n- 



A = 0, (5) 



which generalizes the standard Budden equation by con- 
sidering two distinct resonance layers (at x = ± a) in- 
stead of a single resonance. A similar extension of the 
standard Budden problem associated with the presence 
of sheared flows in magnetized plasmas leads to three 
resonance layers and a phase-space cavity Q . 
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We solve the extended Buddenproblem (p} by using 
the modular-eikonal analysis [2j, |3J, which involves the 
following single-crossing amplitude transmission and con- 
version coefficients (r a ,/3 a ) = (Td,fld) = (t_,/?_) and 
(T b ,(3 b ) = (t c ,(3 c ) = (r +! /3 + ), where 



exp | 



l^±| 2 ) 



P± = v/1 - t± exp(i<p±) 



(6) 



with the conversion phases <y9± = arg[r(i |t/_(_ | 2 )] and 
\rj ± \ 2 = |?7 ± | 2 /(2fco), where 2k = \{D A , D±}\ denotes 
the ray phase-space Poisson bracket of the dispersion 
functions evaluated at the conversion points. The three- 
wave conversion scenario shown in Fig. [T] exhibits a leaky 
cavity feature in ray phase space that allows energy flux 
to recirculate around the quadrangle abed. To represent 
this process, we introduce the complex-valued recircula- 
tion coefficient (denoted by the heart symbol) 
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expressed as the sum of the propagation phase <5> = <f k dx 
(expressed as the sum of the area 4 koa enclosed by the 
uncoupled quadrangle abed and corrections to the un- 
coupled eikonal phases due to finite coupling constants 
rj± =/= 0) and the net conversion phase 'Avg{j3 a f3 b f3 c f3d) = 
2 (ip + +</?_). Note that the magnitude of the recircula- 
tion coefficient 1^1 increases from to 1 as the coupling 
constants go from to infinity. 

All coefficients (T, C±, R) are calculated by adding the 
amplitudes of each successive completed circuit in Fig. Q] 
and, thus, they exhibit interference effects that depend 
on ^. These interference effects enter through the re- 
circulation coefficient (|7|) in two different ways in Fig. [T] 
First, successive powers of V are added to the zeroth- 
order amplitudes for the transmission T, conversion C+, 
and reflection R coefficients; see the recirculation coef- 
ficient W in the boxed expression in Fig. [T] between re- 
gions a and b. Hence, the net transmission coefficient 
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the net conversion coefficient C+ 
the upward fc-propagating wave is 
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C+ = \r c p b T a (1 + Z> + Z> 2 + •••)! 
= t 2 (1-t 2 )t!|1 9|~ 2 , 

the net reflection coefficient R = l^outlVl^rJ 2 i s 

,2 



R = \r d p c (5 b r a (1 + + 9 2 + •••)! 

= T_(l-T 2 ) 2 |1 - q- 2 . 



(9) 



to 



(10) 



(11) 




Coupling 



Separation Parameter k n CT 



10)1 as a function of fcrjcr and 
n(2)/27r]2. At fixed coupling 



FIG. 2: Conversion coefficient 

\ r l- |/vAo) f° r 7 = 1 and \rj, | = 
constant \fj_\, the conversion coefficient C+ exhibits interfer- 
ence effects as the separation parameter koa is varied and the 
interference phase (|8]) passes through integer multiples of n. 



Second, the recirculation coefficient ([7]) enters into 
the conversion coefficient C_ at lowest order because 
of the interference between the zeroth-order conver- 
sion amplitude — /?* and the one-circuit amplitude 
(T a f3df3cflbe l ' s ')Ta (see the boxed expression below re- 
gion a in Fig. [TJ. Each additional circuit around the 
quadrangle adds successive powers of to the one- 
circuit amplitude. Hence, the net conversion coefficient 
C_ ee |i?° ont | 2 /|A° t | 2 to the downward fc-propagating 
wave is 

C_ = \-P: + (r a p d p c p b e^) T a (l + 9 + <? 2 + ---)| 2 



(1 



V\~ 2 [rl(l 



2 \2 



(1 



[(<? + <?*) -2 P| 2 



(12) 



where we used | - a* + b\ 2 = \a\ 2 + \b\ 2 - (a*b* + ab). 
The modular-eikonal coefficients (|9))-(IT2|) depend on the 
model parameters a and |?7_i_| but not the propagation 
parameter 7, which represents a simplifying property of 
the symmetric-quadrangle model used here. They also 
satisfy the energy-flux conservation law 



1 = T + Cj 



C_ 



R. 



(13) 



which follows from Eq. ^j. 

Figure [2] shows a plot of the conversion coefficient 
(I10p as a function of the separation parameter k^a and 
the coupling constant \rj_ | , for 7 = 1 and \rj + \ = 

[ln(2)/27r]^ (where the standard Budden conversion coef- 
ficient reaches its maximum at j); the other coefficients 
(T, C_,R) exhibit similar interference effects while sat- 
isfying the conservation law (fT3|) . We make three re- 
marks concerning our modular-eikonal coefficients ([5])- 
(fl"2"|) . First, the standard Budden coefficients [l[ are re- 
covered from Eqs. ©-(HU) with r\_ — (so that r_ = 1 
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FIG. 3: Conversion coefficients C± versus separation param- 
eter feocr for 7 = 1 and \rJ + \ — \rj_\ = [ln(2)/27r] ' . The solid 
and dotted lines represent the theoretical formulas (I10|l and 
(1121) . respectively, and the symbols <0 and A represent nu- 
merical results for C+ and C_, respectively. The standard 
Budden conversion Cb corresponding to rj- = is also shown 
(represented by a horizontal set of x symbols). 

and /?_ = = <v>). Second, as |»7_| — > oo, we find C — > 1 
and (T, C+, R) — > as can be seen in Figure [2] (where C+ 
becomes very small as \fj_\ > 1). Third, the use of an 
asymmetric-quadrangle model would only shift the min- 
ima and maxima in the modular-eikonal coefficients and 
preserve the interference effects. 

We note that all interference effects disappear when 
the order of the two counter-propagating waves in Fig. [T] 
is reversed (i.e., we assign negative values to the sep- 
aration parameter k^a) since recirculation is no longer 
possible. In this case, the coefficients (l9))- (|12j) simply be- 
come T = t 2 t 2 , C+ = t\ (1 - r 2 ), R = (1 - rlf, and 
C_ = t? (1 — ) , respectively, which still satisfy the con- 
servation law (fT5|) . This reversed scenario may be viewed 
as a combination of a standard double-conversion Bud- 
den problem, with coefficients T b = T + C_, = C + , 
and Rb = R, and an isolated single-conversion problem. 

We now compare the analytical results presented above 
with the direct numerical integration of the coupled- 
wave differential equations ([3]) with 7 = 1, \J}> | = 
[ln(2)/27r]5 = \rj_\, and the boundary conditions corre- 
sponding to only a left-propagating A-wave with eikonal 
amplitude Aq at x — — 00 (and no other waves). With 



these boundary conditions, the energy-flux conservation 
law (@J yields 

-fco|A„| 2 = k (\A R (+)\ 2 -\A L (+)\ 2 ) + £ 

where Ah > l(+) denote the eikonal amplitudes of the right 
(R) and left (L) propagating A-w&ves at x = +00. 
By comparing this equation with Eq. (113[) . we define 
the numerical coefficients T = |A | 2 /|Al(+)| 2 , C± = 
J±(+)/k \A L (+)\ 2 , and R = \A R {+)\ 2 / \A L (+)\\ The 
conversion coefficients C± are shown in Fig. [3] for 7 = 1 
and \rj + \ = \rj_ \ = ^= (the latter value corresponds to the 

strongest interference effects observed in Fig. [5]), where 
excellent agreement is found between numerical results 
and the analytical conversion coefficients (|10p and (|12p . 
The value of the conversion coefficient Cb correspond- 
ing to the standard Budden problem (with rj^ = 0) is 
also shown and it is clearly seen that the presence of the 
third (counter-propagating) wave reduces the conversion 
to the bulk- ion wave, i.e., C+ < Cb- We note that the 
modular-eikonal analysis 0, Q is valid only when the 
area (approximately 4/coc) enclosed by the quadrangle 
in Fig. Q] is large enough compared to 27r (i.e., the sepa- 
ration a is large compared to the wavelength). 

We conclude that the presence of an energetic-particle 
population may interfere with RF heating and current- 
drive processes in high-temperature magnetized plas- 
mas. The simple one-dimensional model presented here 
for an extended Budden problem involving three cou- 
pled positive-energy waves demonstrates these interfer- 
ence effects and our numerical analysis validates the use 
of the modular-eikonal approach. The co-propagating 
and counter-propagating scenarios associated with a 
negative-energy wave supported by an inverted pop- 
ulation of energetic-particles allow t_ = cxp(+ 7r|r?_ | 2 ) > 
1 and yield negative conversion coefficients C_ , which 
can be expected to give bulk-ion conversion coefficients 
C+ that could exceed unity even as the conservation law 
(fl~3| is preserved. We will investigate this potentially 
new alpha-channeling process, as well as develop more 
realistic models, in future work. 
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